Experimental ionization of atomic hydrogen with few-cycle pulses 
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We present the first experimental data on strong-field ionization of atomic hydrogen by few-cycle 
laser pulses. We obtain quantitative agreement at the 10% level between the data and an ab initio 
simulation over a wide range of laser intensities and electron energies. 



The interaction of intense few-cycle infrared laser 
pulses with matter induces tunneling ionization and sub- 
sequent quantum dynamics of freed electrons. Intense 
few-cycle pulses are difficult to generate and use because 
of the stringent requirements on dispersion control over a 
broad bandwidth. However, they offer unparallelled op- 
portunities to reveal and control the electronic dynam- 
ics of atoms [U [5] and molecules [3J 0] and to gener- 
ate isolated attosecond pulses in the extreme ultraviolet 
[5]. The few-cycle regime is particularly challenging for 
simulations, as intensities approaching 10 15 W/cm 2 can 
be reached before the ionisation response saturates. At 
these intensities, a photoelectron driven by intense long- 
wavelength radiation can travel a distance hundreds of 
times larger than the size of the parent atom and can 
have energies of many tens of eV, imposing stringent re- 
quirements on the simulation grid. Ab initio simulations 
in this regime can be carried out only for atomic H due 
to its simple electronic structure. 

Here we describe an experiment on the interaction of 
intense few-cycle laser pulses with atomic hydrogen (H), 
the simplest of all atomic systems and the traditional 
test case for atomic physics. No data on H has previ- 
ously been available in this regime of laser interaction. 
Previous strong-field experiments with atomic H [51 [7] 
used relatively short-wavelength pulses that were many 
optical cycles in duration with maximum intensities of 
10 14 W/cm 2 . Our data show excellent quantitative agree- 
ment, at the 10% level, with ab initio simulation over a 
wide range of electron energies and laser intensities. 

The experimental apparatus is composed of an atomic 
H beam interacting with a few-cycle strong-field laser 
(Fig. 1). The laser used is a commercial Femtolaser 
'Femtopower Compact Pro'. Each pulse has energy 
of 150 /iJ and the pulse repetition rate is 1 kHz. The 
spectral width of the laser is 150 nm at full width half 
maximum (FWHM) centered at 750 nm. The pulse du- 
ration at FWHM of the intensity envelope is 6.3 ± 0.2 
fs at the interaction region, or alternatively ~ 2.5 opti- 



cal cycles. An off-axis parabolic mirror of 750 mm focal 
length is used to focus the beam to a spot size of 47 /im 
1/e 2 radius. The laser carrier-envelope phase was not 
stabilized in these experiments. 

The atomic H beam is created via collisional dissoci- 
ation in a radio frequency (RF) discharge powered by 
a helical resonator [8]. An RF signal at a frequency of 
75 MHz and power of 8 W is applied to the resonator 
and the dissociation efficiency is determined via emis- 
sion spectroscopy of the discharge. The atomic beam 
emerging from the discharge is 80 ± 15% H atoms by 
number with the remainder being undissociated H2 . The 
atomic H beam passes through two apertures, produc- 
ing a uniform-density beam of 0.5 mm diameter and 
< 2 x 10 -3 mrad divergence angle at the interaction re- 
gion. 

The electron detection system is composed of a series of 
electrostatic lenses that act to repel low energy electrons. 
Electrons passing the repeller are accelerated to ~ 250 eV 
before being detected by a channeltron. The channeltron 
is positioned in line with the laser polarization direction 
and is operated in electron counting mode. 

The experimental procedure involves measuring the 
electron yield over 10,000 laser pulses for a range of laser 
intensities and repeller voltages. To obtain the electron 
yield for atomic H alone (Ch), we take electron counts 
for three different settings: with the dissociator on (C on ), 
with the dissociator off (C ff) and a background measure- 
ment taken with the atomic beam blocked (Cbkg)- The 
parameters are combined according to 

Cr = (C on — Cbkg) — (1 — M)(C ff — Cbkg), (1) 

where y, is the fraction of atomic H in the beam. The 
procedure is repeated for a range of peak intensities from 
1.2 to 5.4 x 10 14 W/cm 2 while varying the repeller voltage 
from 10 V to a maximum of 80 V. The peak intensity is 
changed by inserting pellicle beamsplitters, which have a 
negligible effect on pulse duration. The laser intensity is 
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FIG. 1. A beam of H atoms interacts with a few-cycle laser 
beam and photoelectrons are ejected with a wide range of 
energies. Electrons with energy above a cutoff value emerge 
from the repeller and are detected by a channeltron. 



independently calculated prior to each data set by mea- 
suring the optical power using a thermal power meter, the 
pulse duration using an auto correlator and the focused 
spot size using a charged-coupled device (CCD) beam 
profiler. The error in the independent measurement of 
absolute peak intensity is ~ 10%, which is comparable to 
the state of the art [9]- 

As seen from Eq. ([I]), quantitative characterization 
of the dissociation efficiency /i is needed in order to iso- 
late the electron yield due to hydrogen. We use emis- 
sion spectroscopy of the discharge to obtain an accurate 
value for fi. The relative intensities of two atomic H 
lines, Balmer-a and Balmer-/?, and one molecular line, 
(2-2) Ql of the H 2 Fulcher-a system, give the degree of 
dissociation [TU]. Typical values of the intensity ratios 
lB a /I(2-2)Q\ and IbJI{2-2)qi are 393 ± 39 and 74 ± 7 
respectively. A measurement of the rotational series of 
the (2-2) Q Fulcher-ev system must also be made in or- 
der to calculate the effective gas temperature inside the 
discharge [TT]. The gas temperature inferred from those 
measurements is 480 ± 50 K. All intensities are mea- 
sured using aim double monochromator with 0.15 nm 
resolution with detection by a photomultiplier tube. The 
measured dissociation fraction of [i = 80 ± 15% atomic 
H by number is comparable to the dissociation fraction 
measurements in the literature [121 113] . The estimated 
error in the dissociation fraction measurement is inferred 
from errors in the relative intensities of the spectral lines. 

Our experimental data are compared to theoreti- 
cal predictions obtained by direct integration of the 
non-relativistic time-dependent Schrodinger equation 
(TDSE). Two different TDSE calculations were imple- 
mented independently. Both used the velocity form of 
the interaction Hamiltonian H and expanded the wave- 
function in spherical harmonics "J = Y](1" Rg(r)Ya)(0), 
with radial functions Ri{r) defined in a box of size i? max . 
The first method used the Lanczos propagator [14 . This 
method represents the wavefunction as a linear combi- 
nation of vectors "J, H^Sf, . . . , H m ^f , where m = 5. Re- 
sults converged well with values of £ max = 25 a.u. and 
R max = 2000 a.u. and time-step A = 0.01 a.u. The 



second approach was the 'matrix iteration method' |15j . 
Its specific application to the present problem has been 
described recently [TB]. Results converged for £ max = 30 
a.u., i? max = 3000 a.u., and A = 0.04 a.u. Comparison of 
the predictions from the two TDSE approaches resulted 
in excellent agreement and only one set of TDSE simu- 
lations is used in the following. 

The theoretical simulations give momentum-resolved 
electron probability distributions for a fixed laser peak 
intensity. To simulate the experimental signal, we in- 
tegrated the theoretical distributions over the detector 
acceptance function, the Gaussian laser intensity distri- 
bution, and the atomic density distribution in the in- 
teraction region. Careful modeling of the detector ac- 
ceptance and the laser intensity distribution was critical 
to obtain good agreement between experiment and the- 
ory; for instance, approximating the laser intensity by 
a top-hat profile led to disagreement at the 50% level. 
The detector acceptance function was found by simulat- 
ing the electron trajectories in the electrostatic fields of 
the detection apparatus for a variety of trajectory initial 
conditions and electrode voltages. The acceptance angle 
was f=a 7° with slight variations as a function of repeller 
voltage Vd- The electron cutoff energy Ed was found to 
follow the approximate formula Ed w 0.8(Vd — 5), where 
Ed has units of eV and Vd has units of V. 

To account for the Gaussian intensity variation of the 
laser beam, the atomic density is taken to be constant 
within an infinitely long cylinder oriented perpendicular 
to the laser beam direction and zero outside this cylinder. 
The Rayleigh range of the focused laser beam is much 
larger than the atomic beam width, which is again much 
greater than the focused spot size. The overall electron 
yield is therefore equal to the integral of the detection 
probability over the transverse Gaussian profile of the 
laser beam [SJ. 

Figure 2 shows the detected ionisation signal from 
atomic H compared to the TDSE prediction. The atomic 
H electron yield is plotted as a function of electron cut- 
off energy Ed for a number of laser intensities, using the 
approximate linear relationship between repeller voltage 
Vd and cutoff energy Ed- The theoretical predictions 
are compared to experimental data by a fitting routine 
that takes each independently measured intensity mea- 
surement as an input. The full set of data points is simul- 
taneously fitted using only two fit parameters: an overall 
scaling factor and an intensity scaling factor. The con- 
stant overall scale factor applied to the yield accounts for 
the absolute detector efficiency and the target density in 
our apparatus, which are not independently measured. 
The intensity scaling factor gives a fit value for the abso- 
lute peak intensity, while leaving the relative intensities 
of the various data runs at their independently measured 
values. 

As shown in Fig. 2, we obtain agreement at the 10% 
level between the experimental data and the ab initio 
TDSE prediction over a wide range of electron energies 
and laser intensities. Our analysis of experimental un- 
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FIG. 2. Experimental data (dots) versus TDSE predic- 
tions (solid line). The laser intensity range is 1.2 to 5.4 
xlO 14 W/cm 2 and is indicated below the corresponding data 
run. The difference between two TDSE predictions was much 
smaller than our experimental error. 

certainty includes contributions from the dissociation ef- 
ficiency measurement, the shot noise, and an Allan de- 
viation analysis of the electron yield. The uncertainty 
in determining fi results in a contribution of ~ 5%, the 
shot noise is 4 - 10% depending on yield, and an Al- 
lan deviation analysis shows that an additional ~ 5% 
change in electron yield can be expected over the dura- 
tion of data acquisition. The data run with an intensity 
of 3.2 x 10 14 W/cm 2 was a slight outlier and a correc- 



tion factor was found by fitting the run independently. 
The measured intensity was therefore adjusted by 6%, an 
amount consistent with experimental error in the inde- 
pendent measurement of intensity. The excellent agree- 
ment between data and predictions over a wide range 
of parameters provides strong evidence of our ability to 
model the experimental conditions accurately. 

The absolute intensity scale determined by the TDSE 
fit is 0.94 ± 0.01 times the independently measured ab- 
solute intensity scale. Since this scale factor is consistent 
with 1 to within our absolute intensity calibration error 
of 10%, we see that the TDSE fit is capable of extract- 
ing the true peak intensity from the data at this level 
of uncertainty. The small error in the intensity scaling 
for the TDSE fit offers the possibility of laser intensity 
calibration at the 1% level in future experiments. 

We have presented experimental results on few-cycle, 
strong-field interaction with atomic H, the classic test 
system for atomic physics. The excellent quantitative 
agreement with the TDSE model indicates the accuracy 
of our model for our experimental conditions. Similar 
experiments can be accurately calibrated using atomic 
hydrogen and subsequently used with other atomic and 
molecular species. Future measurements of fully re- 
solved electron momentum distributions and HHG spec- 
tra for atomic H will provide benchmarks for strong-field 
physics. 
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